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■ Abstract. A major paradigm shift has recently revolutionized our picture of glob- 
• ' ular clusters (GC) that were long thought to be simple systems of coeval stars born 

i-J^ ' out of homogeneous material. Indeed, detailed abundance studies of GC long-lived 

low-mass stars performed with 8-lOm class telescopes, together with high-precision 
JL ' photometry of Galactic GCs obtained with HST, have brought compelling clues 

< I ' on the presence of multiple stellar populations in individual GCs. These stellar 

, subgroups can be recognized thanks to their different chemical properties (more 

^ ■ precisely by abundance differences in light elements from carbon to aluminium; 

' ' ' see Bragaglia, this volume) and by the appearance of multimodal sequences in the 

colour-magnitude diagrams (see Piotto, this volume). This has a severe impact on 
, our understanding of the early evolution of GCs, and in particular of the possible 

^ ■ role that massive stars played in shaping the intra-cluster medium (ICM) and 

in inducing secondary star formation. Here we summarize the detailed timeline 
I we have recently proposed for the first 40 Myrs in the lifetime of a typical GC 

, following the general ideas of our so-called "Fast Rotating Massive stars scenario" 

• (FRMS, Decressin et al. 2007b) and taking into account the dynamics of interstel- 

I ' lar bubbles produced by stellar winds and supernovae. More details can be found 

I in Krause et al. (2012, 2013). 
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1. Introduction self-enrichment as depicted from star-to- 

^ ' star abundance variations in light elements: 

Two main scenarios have been proposed in polluters could indeed be either fast rotat- 

the literature to account for multiple stellar ing massive stars (FRMS, e.g. Prantzos & 

populations in GCs. They originally differ Charbonnel 2006; Decressin et al. 2007a,b; 

on the nature of the first generation (IG) Schaerer & Gharbonnel 2011; Krause et 

stars potentially responsible for the GC al. 2013) or AGB stars (e.g. Ventura et 
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al. 2001, 2011), with very different conse- 
quences on the understanding of the dy- 
namical evolution of GCs during the first 
40-100 Myr of their life. Here we summa- 
rize the results presented in Krause et al. 
(2012, 2013) where we revisit the FRMS 
scenario for typical Galactic GCs. 

We consider as "typical" a GC where 
no [Fe/H] abundance variation is observed 
from star to star, like e.g. NGC 6752. 
According to Decressin et al. (2010), the 
initial mass of this GC must have been of 
the order of 9 xlO^ M0, with a half-mass 
radius of ^3 pc, and a star formation effi- 
ciency of 30 %. Such a cluster must there- 
fore have hosted ^ 5700 IG massive stars 
(with masses above 25 M0) under the as- 
sumption of a Salpeter IMF for IG stars 
more massive than 0.8 Mq (note that we 
assume a log-normal IMF for IG and 2G 
stars of lower mass). We take into account 
the dynamics of interstellar bubbles pro- 
duced by stellar winds and supernovae as 
well as the interactions of FRMS with the 
ICM and we describe the kinematics by a 
thin-shell model. We can then put forward 
a description of the time sequence for the 
first 40 Myrs of the evolution of GCs (the 
initial time being set at the assumed co- 
eval birth of IG stars) with three distinct 
phases. 

2. Timeline for the first 40 Myrs in 
the lifetime of a typical GC 

2.1. Wind bubble phase and 2G star 
formation 

During the first few Myrs after the forma- 
tion of the IG, the cluster is strongly im- 
pacted by the fast radiative winds of the 
numerous massive stars. Because their in- 
tegrated energy is much smaller than the 
ICM's binding energy, those winds can not 
lift any noteworthy amount of gas out 
of the GC (see also Krause et al. 2012), 
they will create large hot bubbles around 
massive stars and compress the ICM into 
thin filaments. Additionally in this environ- 
ment "classical" star formation is inhibited 



by the high Lyman- Werner flux (see also 
Conroy & Spergel 2011) during the wind 
bubble phase as well as the subsequent su- 
pernova phase fS I2.2l) . 

However if massive stars have equa- 
torial mass ejections that form decretion 
disks as we originally proposed in the 
FRMS scenario (Decressin et al. 2007b; 
see also Krticka et al. 2011) and as ob- 
served e.g. around Be-type stars (see e.g. 
Rivinius et al. 2001; Towsend et al. 2004; 
Haubois et al. 2012), accretion of pris- 
tine gas may proceed in the shadow of 
the equatorial stellar ejecta. Even if a sub- 
stantial fraction of the power of the dom- 
inantly polewards-directed (because of ro- 
tation) fast wind would be channeled into 
the decretion disks, this would still be in- 
sufficient to unbind them from their parent 
stars. Similarly, close encounters with other 
stars at sufficiently small impact parameter 
are too unlikely to affect the disks signifi- 
cantly on an interesting timescale. Second 
stellar generation (2G) low-mass stars are 
then expected to form due to gravitational 
instability in the discs around IG massive 
stars, those disks being fed both by the H- 
burning ashes ejected by the FRMS and by 
pristine gas. 

The total amount of gas that is made 
available to form 2G stars through this non- 
classical mode and with abundance prop- 
erties similar to those observed today is 
high. Indeed Decressin et al. (2007b) mod- 
els predict that FRMS lose about half of 
their initial mass via equatorial ejections 
loaded with H-burning products during the 
main sequence and the luminous variable 
phase. As a consequence about ~ 1 Myr af- 
ter the birth of the IG, the total mass lost 
through this mechanism by all the FRMS 
is of the order of 10" Mg/Myr. On the 
other hand a similar amount of pristine 
gas is expected to be brought by the ac- 
cretion flow, with mixing proportions de- 
pending on the FRMS orbit-averaged ac- 
cretion rate. 2G stars will therefore form 
with the current light element abundances 
of the self-gravitating disk of their massive 
IG parent star, in agreement with the de- 
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tailed abundance calculations of Decressin 
et al. (2007b) that reproduce nicely the ob- 
served [0/Na] distribution in the typical 
GC NGC 6752. 

The formation of 2G low-mass stars is 
expected to be complete after ^ 9-10 Myrs 
(i.e., roughly the lifetime of a 25 Mq star) 
after the formation of the IG. Note however 
that the duration of the wind bubble phase 
is rather uncertain since it depends on the 
mass limit (and therefore on the lifetime) at 
which stars explode as SNe or turn silently 
into black holes, which value is not settled 
yet (see Decressin et al. 2010 for references 
and further discussion on that issue). 

2.2. Supernova phase 

Krause et al. (2012) showed that gas ex- 
pulsion via supernovae explosions, which 
is generally invoked to drastically change 
the GC potential well and to induce sub- 
stantial loss of IG stars (Decressin et al. 
2007a, 2010; Bekki et al. 2007; D'Ercole et 
al. 2008; Schaerer & Charbonnel 2011; see 
also ij 12.31) . does not work for typical GCs. 
Indeed, while the energy released by SNe 
usually exceeds the binding energy, SNe- 
driven shells turn out to be destroyed by 
the Rayleigh- Taylor instability before they 
reach escape speed; therefore the shell frag- 
ments that contain the gas remain bound 
to the cluster. 

Although the SNe-ejecta will likely be 
mixed with the cold phase of the ICM, 
they are not expected to appear in the 
composition of 2G stars. Indeed once the 
first very energetic SNe explode, the ICM 
becomes highly turbulent and, since the 
Bondi-accretion rate depends on , where 
V is the relative velocity between the star 
and the local ICM, can not be accreted any- 
more onto FRMS decretion disks nor (yet) 
onto any dark remnants that are being 
continuously produced (i.e., stellar black 
holes from stars initially more massive than 
25 Mq stars and later neutron stars from 
initially less massive stars). The SNe phase 
is therefore expected to be sterile as far as 
star formation is concerned. 



2.3. Dark remnant accretion phase 

The FRMS and the AGB scenario strongly 
differ on many aspects of the early dynam- 
ical and chemical evolution of the young 
GCs. They both agree however on the fact 
that GCs must have been much more mas- 
sive initially than today and call for sub- 
stantial loss of IG low- mass stars in order 
to explain todays observed proportions be- 
tween 1 and 2G stars (Carretta et al. 2010; 
Schaerer & Charbonnel 2011). Within the 
framework of the AGB cooling flow sce- 
nario, this powerful event is expected to 
occur before the formation of 2G stars. At 
the opposite, 2G stars are expected to form 
much earlier in the FRMS scenario as de- 
scribed in § 12.11 

Krause et al. (2012) proposed that very 
fast gas expulsion, which is expected to 
unbind a large fraction of IG low-mass 
stars sitting initially in the GC outskirts, 
may be powered by the sudden and con- 
comitant activation of IG dark-remnants 
(i.e., 1.5 Mq neutron stars from stars with 
masses between 10 and 25 Mq, and stel- 
lar mass black holes with 3 Mq from more 
massive stars). Accretion of interstellar gas 
onto these objects can set in only once tur- 
bulence has sufficiently decayed in the ICM 
after the type II SNe have ceased, i.e. about 
35-40 Myrs (corresponding to the lifetime 
of a 9 Mq star) after the coeval birth of IG 
stars. Assuming coherent onset of accretion 
onto all the IG stellar remnants and con- 
sidering that they could contribute to gas 
energy at 20% of their Eddington luminos- 
ity, Krause et al. (2012) showed that a suf- 
ficient amount of energy is released in at 
most 0.06 Myrs (i.e., about half a crossing 
time), which is short enough to avoid gas 
fall back due to Rayleigh- Taylor instabil- 
ity and to eject the ICM together with the 
outer weakly bound IG low-mass stars. 

After the major expulsion event de- 
scribed above, the cluster is expected to 
be devoid of pristine gas. However the 
dark remnants may be reactivated later, 
this time through the eventual accretion 
of slow-AGB winds. When this happens 
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we expect again quick gas expulsion, which 
would inhibit any subsequent star forma- 
tion from the AGB cjccta. This situation 
is alike the possible activation of super- 
massive black holes in nuclear star clus- 
ters and elliptical galax;ies by AGB-winds 
(e.g. Gaibler et al. 2005; Davis et al. 2007; 
Schartmann et al. 2010). 

Importantly, Krause et al. (2012) find a 
limiting initial cloud mass of ~ 10^ Mq for 
the dark remnant scenario to work. Above 
that mass, the cold pristine gas may not 
be ejected, and subsequent star formation 
may proceed in the ICM that would then 
be enriched in SNe cjccta. This can explain 
the fact that the more massive (and atyp- 
ical) GCs (Omega Gen, M22, M54, NGG 
1851, NGC 3201; Carrettaet al. 2010, 2011; 
Johnson & Pilachowski 2010; Marino et al. 
2011; Simmerer et al. 2012, 2013) do ex- 
hibit internal [Fe/H] spread, while less mas- 
sive (typical) ones do not. 



3. Conclusions and perspectives 

After the tremendous paradigm shift that 
has recently revolutionized our view of 
GGs, new observational and theoretical ad- 
vances make the present a particularly ex- 
citing time for the study of the origin, for- 
mation, and evolution of these objects in 
the Milky Way as well as in external galax- 
ies. 

On the theoretical front one of the fun- 
damental issues concerns the physics of 
the interstellar medium and its interactions 
with stars. Also, the formation of 2G stars 
in massive disks fed by both equatorial 
ejections of parent FRMS and accretion of 
pristine gas has to be studied in details. 
Cutting-edge strategies have to be estab- 
lished in order to tackle the open issues 
from an innovative point of view. In partic- 
ular, relevant multi-dimensional MHD sim- 
ulations would greatly help in clarifying de- 
tails of chemical and dynamical feedback of 
IG massive stars and of secondary star for- 
mation. Work is in progress in that direc- 
tion. 
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